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By Clarence A. Brown, Jr. 

SUMMARY 


The longitudinal stability characteristics of a simple infrared 
homing missile have been determined in flight with a rocket -powered model 
by the Langley Pilotless Aircraft Research Division. Static and dynamic 
longitudinal stability derivatives of this cruciform, interdigitated, 
canard-wing missile configuration were determined from an investigation 
using the pulse-rocket technique for a Mach number range of 0.7 to 1.4. 

The average lift-curve slope for the model varied smoothly over the 
Mach number range tested. The large tail length of the model was extremely 
effective in increasing the damping-in-pitch derivative throughout the 
Mach number range tested but the damping, in terms of percent critical 
damping, was approximately 10 to 20 percent because of the large inertia 
of the model. The aerodynamic -center location varied smoothly with Mach 
number with the most forward location occurring near a Mach number of 1.0. 
The maximum shift in the aerodynamic -center location occurred between a 
Mach number of 1.0 and 1.35* This shift was approximately 9*0 inches 
(l.8 body diameters). 


INTRODUCTION 


The accuracy of present-day aircraft rocket armament has been ham- 
pered by launching errors, random dispersions, missile reliability, and 
missile complexity. The accuracy might be improved by incorporating some 
type of homing device that would reduce these errors and dispersions. 

An investigation has been undertaken to develop a simple infrared homing 
system that would be reliable as well as rugged. Such a homing system 
has been developed and is described in references 1 and 2. The basic idea 
of this system involves the use of aerodynamics to help in reducing homing- 
system complications and aid in increasing missile reliability. This was 
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accomplished by using a flicker control to produce roll which operates 
directly from the target position as a primary reference, by using a 
rotating lift vector, and by using the rolling of the missile to scan 
the seeker field of view. 

This paper presents the results from a flight-test investigation, 
using the pulse-rocket technique, to determine the static and dynamic 
longitudinal stability derivatives and the drag data for a missile con- 
figuration similar to the one reported in reference 2. The Mach number 
range covered by this investigation was approximately 0.7 to 1.4. The 
model used in this investigation was flight tested at the Langley Pilotless 
Aircraft Research Station at Wallops Island, Va. 


SYMBOLS 


All coefficients are based on a body diameter of 5 inches (0.4l6 foot) 
and a body cross-sectional area of O.I363 square foot. 


a n/g 

a */ g 

•H/g 


b 

d 


C D 


'^in 


normal-accelerometer reading, g units 

longitudinal-accelerometer reading, g units 

transverse-accelerometer reading, g units 

lD"t 

exponential damping constant, e~ , per second 
body diameter, 0.4l6 ft 
acceleration due to gravity, ft /sec^ 
dynamic pressure, lb/sq ft 


I a 7 a n \ W 

drag coefficient, cos a + — sin a — - — 

V S g I <3Ab 


mi nlrmnn drag coefficient 


lift 


coefficient, ( ^ cos a + — sin aj-^— 
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C N 

C Y 

°R 


M 
P 
t 
R 

V 
W 
a 

a = 


pitching-moment coefficient, 

Pitching moment about center of gravity 
qA b d 

TJ 

normal-force coefficient, — — ! — 

S QA h 


lateral -force coefficient. 


_W_ 

g qA b 


re suit ant -force coefficient corrected for trim, 

011/2 


( C N ' C N trlm ) 2 + ( C * - 


(Note that subscript 


^ " Y trimj 

"trim" denotes trim values of model during rolling.) 

Mach number 

period of oscillation, sec 
time, sec 
Reynolds number 

body cross-sectional area at wing-body juncture, O.1363 sq ft 

velocity of model, ft/sec 
model weight, lb 
angle of attack, deg 


— , radi ans /se c 


57.3 at 

0 pitching velocity, radians/sec 

(f rolling velocity, radians/sec 

dC T 


"a. 


— , per degree 


"“•a 


da ' 


per degree 
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dCm 

■Jd/ 


per radian 


~' m a 


d^m 

2V 


per radian 


MODEL AND APPARTUS 
Model Description 


Sketches of the rocket -powered model used in this test are shown in 
figure 1. Sketches of the control surface and wing surface are shown in 
figure 2. Photographs of the model and model-booster combination are 
shown in figures 3 and 4, respectively. Physical characteristics deter- 
mined by preflight measurements are presented in table I. 

The body of the model consisted of two cylindrical sections 5*0 and 
5.5 inches in diameter, and had a body-length — maximum-diameter ratio 
of 24.58 (fig. 1). The 5 •O-inch-diameter section was a modified 
HPAG rocket motor. The nose consisted of a flat face with a drag-reducing 
conical windshield. The conical windshield was supported by an octapod 
mounted in front of the flat nose. Protruding in front of the conical 
windshield was a sting used to mount the angle -of -attack indicator 
(figs. 1 and 3). The canard surfaces were mounted on the 5*5 >-inch - 
diameter cylindrical section of the body, were of tapered plan form, and 
had a maximum thickness at the wing-body juncture of 5-1 percent of the 
chord (fig. 2). The canard surfaces in the horizontal plane were fixed 
at an angle of incidence of 3*95° and the canard surfaces in the vertical 
plane were fixed at zero angle of incidence. 

The 60° delta cruciform wings were mounted on the 5-inch-diameter 
cylindrical section of the body and were interdigitated 45° to the canard 
surfaces. The wing had a modified hexagonal airfoil section with a con- 
stant thickness corresponding to a thickness of 1.4 percent of the chord 
at the wing -body juncture. 

The model used in this investigation differed from the model of 
reference 2 as follows: 

(l) An angle -of -attack indicator was added ahead of the conical 
windshield . 
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(2) The front end of the model was fixed to the rear section of the 
model in such a manner as to eliminate the rolling of the front section 
with respect to the rear section. 

(3) The 5 ^-inch-diameter section of the fuselage was lengthened by 
10 inches but the model overall length was not changed. 

(4) The tracking flares placed on two of the rear fins were 
eliminated. 


Instrumentation 


The model was equipped with an NACA six-channel telemeter which 
transmitted a continuous record of normal, transverse, and longitudinal 
accelerations; angle of attack; total pressure; and rate of roll. The 
transverse, longitudinal, and normal accelerometers were located so as 
to be near the center of gravity of the model when the sustainer motor 
had burned out. Angle of attack was measured by a free-floating vane 
mounted on a sting which protruded ahead of the drag-reducing conical 
windshield (figs. 1 and 3)* Total pressure was obtained by a total- 
pressure tube extended from the fuselage ahead of the wings and in a 

1 ° 

plane 22— to the main wing and canard surfaces. The rate-of-roll instru- 


ment was located just ahead of the center of gravity of the model when 
the sustainer had burned out. 


Model velocity was obtained from the CW Doppler velocimeter and the 
model trajectory was determined through use of a NACA modified SCR 5&4 
radar tracking unit. A radiosonde, released at the time of flight, was 
used to obtain atmospheric data throughout the altitude range traversed 
by the model. 


Test Technique 

The model was launched at 58° 17' elevation angle from a zero-length 
launcher as shown in figure 4. The model was boosted to a Mach number 
of 0.7 by a modified HVAR rocket motor which delivered approximately 
7,000 pounds of thrust for 1.0 second. After separation from the booster, 
a sustainer motor, made as an integral part of the model, delivered 
approximately 2,500 pounds of thrust for 2.6 seconds and propelled the 
model to the peak Mach number of 1.54. After the sustainer burnout, the 
model was disturbed in pitch by a series of six small rocket motors pro- 
viding thrust normal to the longitudinal axis of the model and located 
near the nose of the model. These rocket motors were timed to fire in 
sequence during the decelerating portion of the flight. Transient 
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responses of the resulting disturbances of the model were continuously 
recorded in the form of time histories as the model decelerated through 
the Mach number range. The methods for obtaining the longitudinal sta- 
bility data from the transient responses are presented in references 3 
and k. 


PRECISION OF DATA 
Correction 


The velocity data, as obtained by the CW Doppler velocimeter, were 
corrected for flight-path curvature and wind effects at altitude. The 
magnitudes and direction of these winds were determined by tracking the 
radiosonde balloon. 

In order to obtain the angle of attack at the center of gravity of 
the model, the angle of attack measured at the nose was corrected for 
model pitching velocity and flight-path curvature by the method presented 
in reference 5 - 


Accuracy 

It is believed that the absolute accuracy of the quantities, based 
on the accuracy of the model and ground-instrumentation calibration, are 
within the values tabulated in the following table for two Mach numbers: 
(The magnitude of the random error can be seen by the scatter of the 
points on the curves.) 


M 

Limit of accuracy of - 

M 

cl, deg 

C L 


°Y 

0.80 

0.01 

0.50 

0.26 

0.15 

0.10 

1.35 

.01 

• 50 

.12 

• 07 

.04 


*These values may be positive or negative depending on the 
model and ground-instrumentation zero calibration. 


Parameters dependent upon differences in measured quantities or 
slopes such as Cj_^ are much more accurately determined than the pre- 
viously mentioned errors would indicate. 
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RESULTS AND DISCUSSION 


Data were received for the model tested for a Mach number range 
of 0.75 to 1.40. The Reynolds number of this test ranged from approxi- 
mately 4 x 10^ to 9 x 10 ^, per foot. Variation of Reynolds number with 
Mach number for the test is shown in figure 5* 

Time histories of the lift coefficient Cl, lateral-force coeffi- 
cient Cy, rolling velocity 0, and Mach number are presented in figure 6 . 
Although the model had no controls to produce rolling of the model or a 
roll-control system to prevent rolling, figure 6 shows that as the sus- 
tainer motor burned out (4.7 seconds), the roll rate of the model exceeded 
the instrumentation limit of 10 radians per second. As the Mach number 
decreased (increased time), the roll rate became less and at a Mach num- 
ber of 0.75 the roll rate of the model was near zero. As may be noted 
in figure 6 , irregularities in C-q and Cy occurred during the flight 

of the model. These irregularities of C^ and Cy are due to the 

rolling of the model since the disturbing force occurred only in the 
pitch plane. The irregularities of the lift coefficient made it impos- 
sible to obtain the model period or the exponential damping constant by 
the direct analysis of the time histories of the lift coefficient. 

In order to obtain the period and the exponential damping constant 
of the model, it was necessary to analyze the time histories of the 
resultant -force coefficient by the method presented in reference 4. This 
method consisted of plotting Cjj against Cy for each of the pulse - 

rocket disturbances and, after accounting for the trim as well as pos- 
sible, developing time histories of Cp. A typical plot of the time 

history of resultant -force coefficient is presented in figure 7 * 


Lift Coefficient 

Plots of lift coefficient against angle of attack are shown in fig- 
ure 8 for various Mach numbers. The lift coefficient against angle of 
attack showed some nonlinearities; however, not enough data were avail- 
able to determine any consistent variation with angle of attack. Since 
angle of sideslip was not measured on this model, it was not possible to 
make plots of Cy against 3 or Cp against resultant angle. 

Presented in figure 9 is the variation of average lift -curve slope 
with Mach number. The lift -curve slope was linear within the limits of 
the data and varied smoothly with Mach number over the entire Mach number 
range tested. 
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Dynamic Stability 

The exponential damping constant b is presented in figure 10 for 
the model tested. The rolling of the model and the interaction of the 
normal and transverse motions made it necessary to obtain the exponen- 
tial damping constant from the time histories of resultant -force coef- 
ficient rather than from the angle -of -attack or normal-acceleration 
traces. The method of obtaining the exponential damping constant is 
presented in reference k. Some of the scatter in b can probably be 
attributed to the irregularities in trim which were caused by the rolling 
motions . 

The damping-in-pitch derivative C m( ^ + was obtained from the 

faired values of b (fig. 10) and the faired values of (fig- 9) 

and is presented as a function of Mach number in figure 11. Unpublished 
data for a similar high-fineness-ratio missile configuration show that 
the damping-in-pitch derivative C m( ^ + was of the same order of 

magnitude as that of the present test. The value of for the 

present test is about 10 times that for the model of reference k which 
had about the same ratio of body cross-sectional area to wing area. 

For the Mach number range tested, the large tail length was extremely 
effective in increasing the damping-in-pitch derivative Omq + 0^ but 

the damping, in terms of percent critical damping, was approximately 10 
to 20 percent because of the large inertia of the model. 


Static Stability 

The longitudinal period of oscillation of the model obtained by using 
the time histories of the resultant -force coefficient is presented in fig- 
ure 12 as a function of Mach number. The pitching -moment derivative Cm a 

was obtained from the faired values of the period of oscillation of the 
model ^assuming linear with and is presented in figure 13 as a 

function of Mach number. The pitching-moment derivative (fig. 13) was 
nearly constant throughout the test Mach number range. Preliminary esti- 
mates indicated that the pitching moment would be greater at supersonic 
speeds than at subsonic speeds. Some of this difference between estimates 
and flight data may be attributed to the influence of the flat nose on the 
body of the model. 

Aerodynamic -center location was determined from the curve and 

the faired values of and is presented in figure 14 in terms of inches 

from station 0 against Mach number. Also included in figure 14 are the 
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center-of -gravity locations of the model with sustainer loaded and empty. 
The aerodynamic center moved toward the nose of the model as the Mach 
number increased from 0.75 "to 1.0 and then started a gradual rearward 
movement with increased Mach number. This rearward movement resulted in 
a maximum shift of 9 inches (l.8 body diameters) between a Mach number 
of 1.0 and a Mach number of 1.35* The aerodynamic-center location did 
not appear to have been affected by the rolling of the model. 


Drag 

Drag data for the model tested is presented in the form of C n . , 

•4nin 

based on the fuselage cross-sectional area, against Mach number in fig- 
ure 15 . The drag polars for this model were obtained while the trans- 
verse acceleration was near zeroj therefore, no appreciable drag occurred 
due to angle of sideslip. As may be seen in figure 15, the values 
of ^D min are l ar 6 e ; tut for a configuration such as this these values 

appear to be reasonable. The configuration used in this test was not 
an optimum configuration for drag. As stated in reference 2 this con- 
figuration was chosen in order to use standard components, to have sim- 
plicity in operation, and to keep development cost and tests to a minimum. 


CONCLUSIONS 


The results of a flight investigation of a simple infrared homing 
missile conf iguration for a Mach number range of 0.7 to 1.4 indicated the 
following conclusions: 

1. The average lift-curve slope for the model tested varied smoothly 
with Mach number over the entire Mach number range . 

2. The tests indicates that the large tail length of the model was 
extremely effective in increasing the damping-in-pitch derivative through- 
out the Mach number range tested but the damping, in terms of percent 
critical damping, was approximately 10 to 20 percent because of the large 
inertia of the model. 


10 


NACA RM L 56 DIO 


3 . The aerodynamic-center location varied smoothly with Mach number 
with the most forward location occurring near a Mach number of 1.0. The 
maximum shift in the aerodynamic -center location occurred between a Mach 
number of 1.0 and 1.35* This shift was approximately 9-0 inches (l .8 body 
diameters) . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 5 , 1956. 
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TABLE I 

PHYSICAL CHARACTERISTICS DETERMINED BY PREFLIGHT MEASUREMENTS 

Wing: 

Total wing area in one plane including body- 

intercept, sq ft 2.00 

Wing mean aerodynamic chord, ft 1.24 

Thickness- chord ratio at body juncture 0.0l4 

Wing span, ft 2.15 

Leading-edge sweep, deg 60 

Canard control surfaces: 

Exposed canard area in one plane, sq ft 0.3H 

Exposed canard mean aerodynamic chord, ft 0.481 

Thickness -chord ratio at body juncture 0.051 

Control- surface span, ft 1.125 

General: 

Maximum body diameter, in 5*5 

Body diameter at wing -body juncture, in 5-0 

Fineness ratio 24.58 

Maximum body cross-sectional area, sq ft O .165 

Body cross-sectional area at wing -body juncture, sq ft . . . O.I 363 

Weight, model sustainer loaded, lb 158.5 

Weight, model sustainer empty, lb 116.5 

Moment of inertia - 

X Y> model sustainer empty, slug-ft^ 44.72 

I-£, model sustainer empty, slug -ft ^ 0.20 

Center-of -gravity location, model sustainer empty, 

in. from nose 69.99 

Center-of -gravity location, model sustainer loaded, 

in. from nose 76.88 

Ratio of span of control surfaces to span of wings 0.52 



Model rotated 45° 


Figure 1.- Sketches of model tested. All dimensions are in inches. 
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Control surface 



Wing 

Figure 2.- Sketches of control and wing surfaces for model tested. 

dimensions are in inches. 


All 


Side view 


L-884.91.1 



Model rotated 45° L—884.90#l 

Figure 3.- Photographs of model tested. 
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Figure 4.- Photograph of model and booster prior to launching. 
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Figure 5*- Variation of Reynolds number, per foot, with Mach number . 
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Figure 6.- Time histories of model tested. 
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Figure 6.- Continued 
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Figure 6.- Continued. 
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Figure 6.- Concluded. 
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Figure 7*- Sample time history of re suit ant- force coefficient. 
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Figure 8.- Variation of lift coefficient with angle of attack. 
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Figure 9.- Variation of the average lift-curve slope with Mach number 
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Figure 10.- Variation of the exponential damping constant b 

number. 



with Mach 


Figure 11.- Variation of the aerodynamic damping- in-pitch derivative 

C™ + with Mach number. 
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Figure 12.- Variation of period of oscillation with Mach number. 
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Figure 13.- Variation of the static stability derivative Cn^ with Mach 


number 


Model station in inches 
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Figure 14.- Variation of the aerodynamic- center location with Mach number 
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Figure 15 •- 


Variation of the minimum drag coefficient with Mach number . 
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